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We report a detailed investigation ol the interplay between size quantization and local scattering 
centers in graphene nanoribbons, as seen in the local density of states. The spectral signatures, 
obtained after Fourier transformation of the local density of states, include characteristic peaks 
that can be related to the transverse modes of the nanoribbon. In armchair ribbons, the Fourier 
transformed density of states of one of the two inequivalent sublattices takes a form similar to that 
of a quantum channel in a two-dimensional electron gas, modified according to the differences in 
bandstructure. After addition of the second sublattice contribution, a characteristic modulation 
of the pattern due to superposition is obtained, similar to what has been obtained in spectra due 
to single impurity scattering in large-area graphene. We present analytic results for the electron 
propagator in armchair nanoribbons in the Dirac approximation, including a single scattering center 
within a T-matrix formulation. For comparison, we have extended the investigation with numerics 
obtained with an atomistic recursive Green's function approach. The spectral signatures of the 
atomistic approach include the effects of trigonal warping. The impurity induced oscillations in the 
local density of states are not decaying at large distance in few-mode nanoribbons. 

PACS numbers: 73.22.-f,73.22.Pr 
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I. INTRODUCTION 

In graphene, scattering centers such as impurities, defects, adatoms, and substrate inhomogeneities greatly influence 
the local electronic propertiesii In some samples, the material quality is so high that a single or a few such scattering 
centers can influence the whole device. This may degrade device function, but can also be taken advantage of by 
making various sensing devices?^ Great attention has therefore been focused on understanding the influence of 
scattering on the electronic properties of graphene.— 

In this context, the scanning tunneling microscope (STM) is becoming of increasing importancei^ By utilizing its 
various modes of operation, the STM can be used to map out topography, local density of states, local charge density, 
and more. In this way, a variety of properties of graphene have been revealed. A few examples include perturbations 
in the local density of states around impuritiea^i'^ or near step-edges in the substrate;^ charge puddle formation caused 
by molecules trapped between graphene flakes and the Si02 substrate)^ and resistance caused by step o^^'^^ in the 
substrate or multilayer regions^^ in epitaxial graphene on silicon-carbide. 

At the same time, encouraging progress has been achieved with fabrication of graphene nanostructures. Top-down 
approaches include nano- lithography^ scanning probe methodsfi^ etching with metal nanoparticles along certain crys- 
tal directionsji^ and utilization of the transmission electron microscope (TEM) to simultaneously image and sculpture 
graphene.— A bottom up approach based on chemical synthesis has also been demonstrated.— Another approach in- 
volves unzipping of carbon nanotubesJ^ With that method, the theoretically predicted zero-energy (midgap) edge 
states of nanoribbons with zigzag edge a^°i^^ were directly mapped out by scanning tunneling spectroscopy (STS)<^ 
Theory also predicts that by controlling the width and edges of nanoribbons, a bandgap can be opened up at the Dirac 
point through quantum confinement (see the review [23[). With further progress it may soon become possible to study 
in much greater detail the interplay between quantum confinement and impurity scattering in graphene nanoribbons. 

Many theoretical studies of graphene nanoribbons have been reported in the literature, see the collection of review 
articles in Ref. [H. The effect of impurity scattering and the effects of edge disorder on electron transport have been 
reported in several numerical works. In an effort to simulate the typical experimental situation, random disorder 
is included and the scaling behavior of resistivity with length of the ribbon is studied, revealing different transport 
regimes depending on ribbon width and disorder properties. Here, we go back to the well defined problem of a single 
impurity in order to study in detail the effects on the FT-LDOS. 

In this paper we present results for the spectral signatures of a local scattering center in graphene, taking into account 
quantum confinement in a nanoribbon geometry. This study generalizes the consideration of FT-LDOS of a single 
impurity in bulk graphen o^^i^^ to the case of nanoribbons. We focus the analytic analysis on armchair ribbons in the 
Dirac approximation (linearization around the K-points in the graphene bandstructure), for which the wavefunctions 
and propagators for clean ribbons are known, and solve the impurity problem in a T-matrix formulation. Thereby 
we obtain the electron propagator for an armchair nanoribbon including the effects of a local scattering center. The 
Fourier transformed density of states (FT-LDOS) is then obtained and explained in terms of scattering processes of 
Dirac quasiparticles confined in the ribbon. We extend the analysis to an atomistic tight-binding model of graphene, 
utilizing a numerical recursive Green's function approach. The main effect of going beyond the Dirac approximation 
is trigonal warping, which shows up as a triangular distorsion of the FT-LDOS patterns. 

For comparison we include an analysis of the FT-LDOS in a quantum ribbon in a two-dimensional electron gas 
(2DEG). Many features of the FT-LDOS patterns in graphene ribbons can be understood from the somewhat simpler 
case of a 2DEG, and the new features special for graphene can be highlighted. These include a more complicated 
bandstructure due to the two inequivalent K-points, trigonal warping, as well as interference effects due to the bipartite 
lattice of graphene. 

The outline of the paper is the following. In Section |ll] we discuss the Fourier transform scanning tunneling 
spectroscopy method and illustrate the basic scattering processes at play in a nanoribbon. In Section IIIII we present 
results for the FT-LDOS in a 2DEG quantum channel. In Section ITVl we report our results for the FT-LDOS in an 
armchair graphene nanoribbon within the Dirac approximation and compare with the 2DEG case. In Section |V] we 
present results of numerical simulations of a tight-binding model, including also zigzag nanoribbons as well as effects 
of edge disorder on the FT-LDOS. In Section IVTl we summarize the paper and give some conclusions and an outlook. 
Most technical results of the analytic analysis have been collected in the Appendices. 



II. FOURIER TRANSFORM SCANNING TUNNELING SPECTROSCOPY 

A scattering center induces a perturbation of the local density of states in its vicinity. For elastic scattering, the 
impurity scatters electrons between states ki — ^ k2 with eg^ = eg^, i.e. on a contour of constant energy E. This leads 
to interference and a wave pattern in the local density of states near the impurity with wavevectors (f = ^2 — fci . After 
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Fourier transformation of the local density of states p(r, E) J^{q, E), the wave vectors of the interference pattern 
are highlighted. The resulting pattern in J^{q, E) can then be used to infer the band dispersion eg. For instance, this 
has been done for metal surfacesi^ This method has also become a valuable tool for probing the properties of high- Tc 
superconductors 3^ 

It is worth mentioning that, neglecting electron-electron interactions, the interference patterns in the local density 
of states discussed above are related to the Friedel oscillations in the electron density n(r) through integration over 
energy including the Fermi-Dirac distribution function, n{r) = —e J p{r, E) f{E)dE. 

By using the STM, the local density of states can be extracted as function of energy by applying a finite voltage 
between tip and sample, i.e. by employing scanning tunneling spectroscopy (STS). By combining Fourier transfor- 
mation with STS, the band dispersion can be studied in the vicinity of the Fermi energy. This method has therefor 
become a valuable spectroscopic tool sometimes called Fourier transform scanning tunneling spectroscopy (FT-STS). 
In graphene, the Fermi energy itself is tunable by a back gate voltage on the substrate graphene is resting on. 
Thereby, FT-STS is potentially a valuable tool for studies of graphene. Indeed, experiment reproduce the graphene 
bandstructure4 

STS bears similarities with angle-resolved photo-emission (ARPES). STS is ideal for spatially inhomogeneous sys- 
tems, while ARPES relies on large-area spatially homogeneous samples. Indeed, STS measures the spatially resolved 
spectral function, i.e. local density of states p(r;E), while ARPES measures the momentum-space spectral function 
A{k;E). By generalizing STS to FT-STS, i.e. Fourier transforming p{f; E) — )■ Af{q;E), a spectroscopy has been 
introduced that can be used to study materials, although we should remember that Af{q', E) is not equal to A{q; E). 

One advantage of FT-STS is the possibility to study nanoscale systems with high spatial resolution. In this paper 
we will investigate the consequences of quantum confinement on impurity scattering in graphene, as seen in FT-STS. 

In Fig. [T] we display a cartoon of a typical scattering process that contributes to the correction to the local density 
of states in a quantum ribbon with one impurity. For simplicity we here discuss the situation in a 2DEG quantum 
channel. Quasiparticles occupying for instance mode m, propagating in the positive y-direction with wavenumber 
1 passes the probing position r = {x,y), after which they can be backscattered by the impurity at fi into mode 
n with wavenumber k„ and propagate back to the probing position r. In this example we neglect evanescent modes 
for simplicity. The contribution to the full propagator from this scattering event will be proportional to the free 
propagators before and after scattering and the potential strength 7, 

Gnni{r,r;E) = -i—^ ; Xn{x)xn{xi) 

( . M \ e"^"l«-^-l , , , , 

where Xn{x) = \/2/W sin(?T,7ra;/M^) is the transverse wavefunction in mode n > 1, and ^ is the electron effective mass. 
Taking into account multiple scattering by the impurity, the potential strength 7 is replaced by a T-matrix. When 
we take the imaginary part of the propagator to get the local density of states, we get spatially oscillating terms 

CX COs[(k„ + Krn)\y ^ 2/i|], 

and 

CX sin[(K„ + Kyn)\y - yi\], 

since the T-matrix is a complex number due to multiple scattering. After Fourier transformation, we find peaks at 
qy = ±(«;„ -I- Urn) and at qx equal to combinations of n and m times tt /W . Thus, in a FT-STS picture of a quantum 
channel, there will be a discrete number of peaks that reflect the available modes. We also note that the Friedel 
oscillations (neglecting electron-electron interactions) will at low temperature oscillate without decay far from the 
impurity site. 

To probe a 2DEG quantum channel with an STM in the way described here will be challenging since the channel 
is typically hidden deep down in a semiconducting hetcrostructurc. Graphene, on the other hand, is 100 % surface 
and directly accessible. 

III. FT-LDOS: RIBBON IN A 2DEG 

In this Section we improve the above discussion to the general case of multiple scattering in a multimode 2DEG 
quantum channel of width W with a single impurity scattering center at r^. The results of this Section will be 
referenced in the following Sections on graphene in order to highlight the distinguishing features of confined Dirac 
quasiparticles. 
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FIG. 1. An electron quasiparticle initially in mode m, with longitudinal wave vector Urn rnay be backscattered by an impurity 
at fi into mode n. The local density of states at r is changed due to interference of the initial and final waves. This leads to 
an interference pattern around r*;. 

Consider the probability amplitude for an electron in the channel to propagate from one point r' to another point 
r. For free propagation in mode n, the amplitude is given by the free propagator (unperturbed Green's function), 
(7„(f, r '; E). In the presence of the impurity an electron initially in mode m may be scattered into mode n. The effect 
of such an extra process will modify the propagator by adding a second term 

GnAr.r'^E) ^ gn{r,n-E)T{n-,E)grn{n,r';E), (1) 

so that the new Green's function will be 

Gn^{r,f'\E)^g^{f,r',E)5^,^ + Gn^{r,f'-E). (2) 

The factor T{fi]E), see Eq. (|A12p . includes multiple scattering by the impurity. The full probability amplitude for 
propagation from r ' to r is given by summing over all mode indices 

G(f,r";i?) = ^G„™(f,r";i?). (3) 

nm 

We may now proceed with the local density of states (LDOS). The correction to the LDOS by impurity scattering 
can be written as 

p{r\E)^ --y^ Im Gnra {r, r; E) 

TT ^ — ' 

nm 

= --Y^lCr.ra{E)~pl,^{x-E)~pl„^{y-El (4) 

nm 
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where the expressions for the factors JCnm{E), Pnmi^)^ PumiV'-E) are given in Appendix El The Fourier trans- 
formed focal density of states (FT-LDOS) can now be computed as 

M{q; E)=--Y, ICrUEWnrniQ.W^rniQy; E) , (5) 
TT — ' 



where 



°° (6) 



and 



TT W 



KME)= r ^e~^^yypy^{y;E). (7) 

The function p^^{x), originally defined on the interval [0, W], is extended to [— W, W] and assumed to be even with 
respect to the origin. Due to the finite width, of the integration interval, the spectral z-component is fixed to 
be integer multiples of tt/W. This is a trick to be able to resolve the minimum change of transverse momenta, n/W, 
when scattering between two different modes. 

It is important to realize that both propagating and evanescent modes play a role in this scattering problem. The 
longitudinal momentum is , = ^2iJLE/h? ~ {mr/WY, where p is the electron mass and n > 1 is the integer mode 
index. At the bottom of a subband, k„ — >■ 0, and the evanescent mode extends far from the impurity and play an 
important role. On the other hand, for energies far from any subband bottom, the local density of states is only 
affected by the evanescent mode in a small region near the impurity. In the discussion of the FT-LDOS we can then 
safely neglect evanescent modes in the sums in Eq. ([5]). The evanescent modes are still taken into account in the 
scattering processes at the impurity through the T-matrix equation, where intermediate modes can be evanescent, 
while initial and final modes are propagating. In all of our numerical calculations, we include 10 evanescent modes. 
Adding even more evanescent modes does not qualitatively change our results. As have been shown, a delta-shaped 
impurity with a finite number of evanescent modes will model an s-like scattereri^ We can now find the different 
components of the FT-LDOS to be 



and 



where 



2W 



Kmily) = — [S^iHn + Urn + %) + + - Qy)] , (9) 



SV^a) = hm -v^-fh + -^) (10) 



and where <7p/e are positive, independent constants defined in Eq. (|A13I) . 
The factor fCnm (E) is given by 



and depends on the scatterer strength 7 and the transversal wave functions Xn{x) = \fljW sin(7T,7rx/Vl^). 

Together, these components give rise to a number of selection rules that govern the modification of the FT-LDOS 
by impurity scattering. To illustrate, we select a narrow channel (M^ — 50ao, where ao defines the unit length) and a 
low energy (E = 0.2t, where r defines the unit energy), such that only a total of three propagating modes are open. 
The scattering FT-LDOS \M{^q\E')\ for the case of the impurity in the middle of the ribbon (x,; = VF/2,j/i = 0) is 
displayed in Fig.[2Ia). Since the scattering is elastic, energy conservation requires that the transverse and longitudinal 
momenta, both before and after scattering, satisfy the relation 2p.ElY? = + k^{E), which is the circle shown in 
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FIG. 2. (a) FT-LDOS of a 2DEG ribbon of width W = 50ao at energy E = 0.025r, with three open (propagating) modes, (b) 
Energy contour of the 2DEG dispersion. The symbols indicate the aUowed momentum values between which scattering can 
potentially take place (circles, diamonds and squares corresponds to n = 1, n = 2 and n — 3 respectively). The two arrows in 
figure (b) illustrates two possible scattering processes that gives rise to the two encircled dots in figure (a). The impurity is 
placed at Xi = W/2 and the n = 2 subband is not scattered by the impurity because the impurity has been located at a node 
of the corresponding transverse wavefunction. 



Fig.mjb). In the channel, the transverse momentum is quantized, — > kn = mr/W, and the only allowed momentum 
values between which the electrons can scatter are indicated by dots and squares on this circle. The FT-LDOS is 
therefor non-zero only at a few. finite number of (f-points. All of of these points lie inside the dotted circle of radius 
2^/2ilE/W shown in Fig. ^^a). 

The factor K,nm{E) will be non-zero only if the transverse wavefunctions of mode n and m have a finite overlap 
at the position of the impurity. Since we have positioned the impurity at Xi = W/2, ICnm{E) will in this example 
be non-zero only if n and m are both odd integers since all the wavefunctions with even indices will have a node at 
X = Xi. Thus, modes with even number n are not scattered by the impurity in this example. 

To understand the exact locations of the g-points, we start by looking at the case Qx = (i.e., 1 = 0). Since all mode 
indices have to be odd, the term Af^{l = 0) will be non-zero only when n = m, i.e., when (rt = 1,to = 1) or when 
(n = 3,TO = 3). This tells us that the points along qx = arc all due to intraband scattering. The factor J\fy{qy;E) 
peaks when Qy = ±2|ki| or when qy = ±2|k3|. These are the four points we see along the line I = 0. 

When qx = tt/W {I = 1), at least one of the indices n and m will be even, and the factor JCnm{E) is zero. This is 
why we see no bright points along this line. This also happens for I = 3 and 1 = 5. 

Along the line qx = 2Tr/W, we have that J\f'-'{1 = 2) is non-zero only when (n = 1, m = 1), [n = l,m = 3) or when 
(n = 3, TO = 1). The factor M'^{qy] E) peaks at qy = ±2|ki| or when qy = ±|ki + K3I, and we see that we have spots 
at these locations along / = 2 in the figure. 

At ^ ~ 4 we must have [n = 3,to = 1) or {n = l,m = 3), which tells us that qy = ±|k3 + At I = 6, we must 
have (n = 3,TO = 3) and qy = ±2|k3|. A similar argument can be made for I < 0, and we can therefor say exactly 
which scattering processes contribute to each dark spot in Fig.l^fa). 

If the impurity is not located exactly at the middle of the ribbon, the even subbands will also be part of the 
scattering process. This is illustrated in Fig. [31 where we have numbered the subband transitions corresponding to 
each bright point. 

In Fig. m we show the result for a wider ribbon calculated both analytically and by doing a recursive tight-binding 
simulation. The parameters are adjusted such that both cases have 20 propagating modes open, and we sec that the 
main features of our analytical calculation and the numerical simulation coincide. 
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IV. FT-LDOS: RIBBONS OF GRAPHENE 



The procedure of calculating the effect of a single impurity on the local density of states in a graphene armchair 
nanoribbon (AGNR) much follows that used for the 2DEG case. Due to the bipartite structure of the graphene 
honeycomb lattice, the propagator G(r, r ';£■) is a 2x2 matrix in sublatticc space, denoted by A and B. We therefor 
start by finding the impurity contribution to the local density of states on each sublatticc. The resulting expressions 
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FIG. 4. (a) Analytical FT-LDOS of a 2DEG ribbon of width W = 200ao and energy E = 0.05r. (b) FT-LDOS taken from a 
numerical tight-binding simulation of a ribbon of width W = 200ao and energy E = — 3.95r. The energies are adjusted such 
that each ribbon has a total of 20 propagating modes open. 



for the A- and B-sublattice LDOS can be written as 



dc nra ^^2) 



TT 

dc nm 



where G^^J^ (f , r ; E) are the two diagonal components of the propagator matrix Gdncm {r, r; E) . The summation 
over the variables c and d are added to account for scattering between different sets of non-equivalent Dirac cone pairs 
and K^. A further elaboration on this is found in Appendix [Bj together with derivations of the expressions for 

p^^{x;E) and P^^im^{y]E). As discussed in section Hill we only need to sum over propagating incoming and final 
transverse modes, labeled by m and n. 

We compute the FT-LDOS on each sublattice as 

M^'^iT^E) = --Y,Ml„M.-,E)MiT^y{qy-,E), (13) 

TT ^ ^ 



The total FT-LDOS is found as a superposition of the two sublattices 

M{q] E) = M^{q] E) + e-'''"'^y {q; E), (14) 



where the extra phase-shift is introduced since the two sublattices are spatially separated by the carbon-carbon 
distance oq in the y-direction. 

Since the transverse wavefunctions, Xnix) — \/l/W sin(n7r/H^a:;), in our AGNR only differ from those of the 2DEG 
by a factor of l/%/2, we have that M^milx; E) = J\f^^{qx]E)/2, as defined in Eqns. © and ©. 
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FIG. 5. Analytic FT-LDOS, M{q\ E), of an AGNR having width W « 50, S = 0.4|r|, 7 = 10|r| and 60 propagating channels. 



The longitudinal FT-LDOS expressions for each sublattice are given by 

3 



c=l 
2 3 



where Ad„cm(-B) = sgn(£')[Kd„(£;) 



£i=l C=rf+1 

+ S^^'^^y{^,r.dm{E) + + 5(^/^)^(Ad„,™(£;) + qy)] , 

'«dm(£')] + Kd — Ky. The two if- independent constants are found to be 



(15) 




[E] (E) 



(16) 



and IC^ncmiE) = ^(^//l-^l)^^dncm(-^)- ^he transverse and longitudinal momenta are now cone set dependent, 
and changes to kdn = nn/W — Kdx and Kdn{E) = \J [E/vfY — fc^„ respectively. Here, S^y[a) = Sy{a) is the 
same function as used in the 2DEG case and defined in Eq. ([T0| . and S^y{a) = fdncm{E)Sy{a) where fdncm{E) = 
[—knkm + Kn{E)Km{E)] + isgn{E)[k„Km{E) + fc,„K„(i5)]. The S'-terms for the AGNR A-lattice have the exact same 
form as the corresponding terms in the 2DEG, while the B-latticc terms arc scaled by a complex mode dependent 
prefactor. This is a consequence of our choice of impurity potential: an impurity fully localized on one A-atom, see 
Eq. (IbTTI) . 

In Fig. [5l we plot \M{q;E)\ for a semiconducting AGNR of width W ~ 50 nm, where E = OA\t\ such that 60 
channels are propagating. The positions, outer shapes, and sizes of the circular features (of radius 2E/vf) are the 
same as those found when studying a sharp impurity in bulk graphene, and are due to the graphene bandstructure. 
In addition to the bulk graphene features, we also see an added rich inner structure due to transverse confinement 
in the nanoribbon. Each peak corresponds to scattering processes that change the transverse momentum by integer 
multiples of tt/W, and change the longitudinal momentum such that the arguments of at least one of the many 
S'^-terms in Eq. (fT5|) vanishes. 

A schematic illustration of one such scattering process is shown in Fig. [S] An electron, initially in cone pair c = 1 
and mode m, is described by a plain wave with momentum —Kim in the longitudinal direction, and a superposition of 
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FIG. 6. Schematic picture of one possible scattering process in AGNR's. Here, the electron (initially in mode m, represented 
by the red dots), is scattered into mode n (green squares). The FT-LDOS will be finite at the g*- values illustrated by the solid 
{qx = (m ± n)n/W) and dotted {q^ — — (m ± n)TT/W) arrows. 

two plain waves with momenta ±fcim in the transverse direction (see the lower red dots). After scattering (within the 
same cone pair) to mode n, the momenta are changed to ki„ and ±fci„ in the longitudinal and transverse directions 
respectively (see the upper green squares). A Fourier transform of the LDOS is proportional to a product of the 
electron wavefunction before and after the scattering event, where each wavcfunction is a linear combination of two 
transverse parts. The FT-LDOS will therefor be finite at the (f- values corresponding to the four arrows shown in the 
figure. Here, qy = |ki„ + Kiml, and ~ (m ± n)7r/W (solid arrows) or — ~{m ± nj-K /W (dotted arrows). When 
scattering to a different cone pair d ^ c, we instead have qy = |Ac;„im|, sec Eq. (|15p . 

When we zoom in on the circular feature in the middle [shown in Fig. [71[a)] we see that the outer ring of non- 
vanishing g-points in |A^(q';i?)| appears to be attenuated compared with what is seen on e.g. the A-latticc alone 
[|7W"^(q'; £^)| shown in Fig. [Tl^b)]. This is due to destructive interference when adding the A- and B-lattice FT-LDOS 
contribution together, as done in Eq. (jl4p . The g-points on the outer circle comes from scattering processes which 
maximimize the change in momenta while still scattering within the same cone pair, i.e., where d = c and fccm — ^ —kdn 
and vice versa. In this case, we have that fdncm = {E/vf)^ which tells us that Mnmily'i E) = —M^y{qy; E) so that 
when the phase factor e"'"^"" in Eq. ((Ti)) is close to unity, the contributions from the A- and B-lattice will cancel each 
other out. Similar cancellations may be seen in Fig[5l e.g. in the circular features to right and left of the central one. 
For other processes and if-values, the interference between the two lattice contributions may not play an important 
role, or we might have constructive interference instead. 

The FT-LDOS is left-right mirror symmetric around the line qx = 0, see Fig. [S] This symmetry appears because 
for every process adding a component q^ in the FT-LDOS, there is another process adding a component q~ , where 
q~ = —q^, see the solid arrows in Fig.O After summation of all such processes, the FT-LDOS acquires the left-right 
symmetry. 

As a consequence, the feature centered around if = is always mirror symmetric by the above argument . On the 
other hand, there is not necessarily a mirror symmetry within the other circular features (i.e. mirror symmetry with 
respect to the individual cone centers). For metallic AGNRs, the transverse modes are constructed from wavevectors 
symmetrically positioned with respect to the cone center (plus the metallic mode at the cone center). See, for instance, 
the two wavefectors k^n and k^n' = —kdn hi Fig. [9] For semiconducting AGNRs, the wavevectors are not symmetrically 
positioned with respect to the cone center, i.e. kdn' 7^ —kdn for any n' . Therefor, the inner structure of the circular 
features centered at finite q are symmetric for metallic AGNRs and asymmetric for semiconducting AGNRs. This 
is illustrated in Fig. [5] for the semiconducting and metallic cases in (a) and (b) , respectively. We conclude that by 
looking at what symmetries there are in the FT-LDOS, one can extract information about whether or not an AGNR 
is metallic or not. 



V. 



NUMERICAL SIMULATIONS OF THE TIGHT-BINDING MODEL 



For our numerical simulations, we use a tight-binding model described by the Hamiltonian 
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FIG. 7. Zoom-ins of the central circular feature of Fig. [5] showing (a) the combined A- and B-lattice contributions and (b) the 
A-lattice alone. When adding the two lattices together, the outer circle is attenuated by destructive interference between the 
two lattice contributions, M^^^'^^" {qy; E), to the total FT-LDOS. 
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FIG. 8. Zoom-ins of the north-eastern circular feature in Fig. (5] In (a), the AGNR is semiconducting and the left- and right- 
hand side is not mirror-symmetric. In (b), the ribbon is made metallic by removing 4 rows of carbon atoms, which restores the 
left-right symmetry again. 



where cj and Ci arc creation and destruction operators for site i. The onsite energy of site i is denoted e;. and the 
hopping amplitude between sites j and i is denoted . The number of atoms in the system is denoted A^. We assume 
that Tij is always zero except when the sites i and j are nearest neighbours. 
The retarded Green's function matrix is defined as 



(18) 
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FIG. 9. Schematic picture of scattering processes conserving mirror symmetry with respect to = (solid arrows) and with 
respect to = ±Kdx (sohd + dotted arrows). 



where is a small positive number. Even though the Hamiltonian is sparse, when written down as a matrix in 
site index space, direct inversion is not a viable alternative when the number of atoms A'' grows large. Instead of 
direct matrix inversion, we use our own implementation of a recent algorithm^ in which the system atoms are added 
onc-by-onc. in a recursive manner. This allows us to save both memory and time, and once we have found all the 
retarded propagators between the system leads and atom i we can calculate the lesser Green's function, defined as 



where I is the lead number (Z = 1, 2 in the case of a simple ribbon), and a; and /3; are indices running over all atoms 
belonging to the surface of lead I. Here, fi{E) and Y,i{E), are the distrubution function and the self-energy of lead /, 
respectively. 

The local density of states on atom i is found from 

p,,{E) = --ImG<{E), (20) 

TT 

and the FT-LDOS is given by doing a discrete Fourier transform over all system atoms, 

TV 

N 



1 N 

AA(<f;i?) = -^e— >V.(i?), (21) 



where fi is the real space coordinate vector of atom i. 

In Fig. llOf a). the result of such a tight-binding simulation is shown for a ribbon and setup matching the one used in 
Fig. [5l with a delta-like impurity placed in the middle {W « 50 nm, 60 propagating channels and Xi = W/2). Upon 
inspection, we notice that the general features are similar compared with our analytical results. Some points, such as 
the outline of the central circle, are attenuated. The tight-binding ribbon do, however, show clear signs of trigonal 
warping due to the dispersion not being perfectly linear. In Fig. [TOlb). we have moved the impurity to the edge of 
the ribbon and we notice that the resulting FT-LDOS image is not very different from the one with the impurity in 
the middle of the ribbon. In Fig. llOf c). we have made the impurity more gaussian shaped (long-range), which leads 
to suppressed scattering and attenuated features. For bulk graphcne, it is well known that a long range impurity can 
not scatter between valleys. In the FT-LDOS, the features centered at q = K^™ are then absent. This is not the case 
here, since the armchair nanoribbon has only one cone in its band structurci^ 

In Fig. [Til we present results for the FT-LDOS of zigzag graphene nanoribbons (ZGNRs). In this simulation the 
ribbon has N = 468 atoms in its unit cell {W « 50 nm), 7 = lOr and E = 0.45|t|. This gives 35 propagating modes. 
In Fig. Illf a). the impurity is located in the middle of the ribbon and we see a pattern very similar to that of the same 
impurity configuration in an armchair ribbon, but with all features rotated 90 degrees due to the different ribbon 
alignment (for ZGNRs', ky is quantized instead). The result of a single impurity on the edge is shown in Fig. [TlTb). 
and in Fig. lllf c) we show the spectra for a ribbon also having rough edges. In the last figure, Fig.[TT](d), we have used 
a gaussian shaped (long-range) impurity, and we here see clearly that inter-valley scattering is now fully supresscd. 
Indeed, since the ZGNR has two cones in its bandstructure, this case is similar to bulk graphene. 



VI. SUMMARY 



In summary, we have presented results for the FT-LDOS of graphene nanoribbons with local scattering centers. The 
interplay between size quantization and scattering leads to characteristic peaks that can be related to the transverse 



13 




03 
D" 



c) 



-3 

■a 

-2 '• 
-1 



1 - 

2- ^ 
3 

4: ^ 



{ ) 



-400 



b) 



-400 -200 200 

q W/ji 



400 



-200 200 400 

q W/ji 

'y 



FIG. 10. Numerical tight-binding FT-LDOS of three AGNRs' {N = 810 atoms in the unit ceU) with diflterent impurity 
configurations, (a) Single impurity, (b) edge impurity, (c) smooth (Gaussian long-range) impurity, where ~ 50 nm, 7 = 10|r|, 
E = 0.4|r| and 60 propagating channels. 



modes of the nanoribbon. The main features include ring-like structures, analogous to the case of an infinite 2D 
graphene sheet with a single scattering center. Inside the ring-like structure, new peaks appear that are related to 
inter and intra band scattering in the ribbon. We have presented analytic results for the electron propagator in 
armchair nanoribbons in the Dirac approximation, including a single scattering center within a T-matrix formulation. 
We have also extended the investigation with numerics obtained with an atomistic recursive Green's function approach. 
The spectral signatures of the atomistic approach include the lifting of degeneracies of transverse modes in the Dirac 
approximation, as well as effects of trigonal warping. The impurity induced oscillations in the local density of states 
are not decaying at large distance in few-mode nanoribbons. 
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FIG. 11. Numerical tight-binding FT-LDOS of four ZGNRs' (iV = 468 atoms in the unit ceU) with different impurity config- 
urations, (a) Single impurity, (b) edge impurity, (c) rough edges, (d) smooth impurity (Gaussian long-range), where ~ 50 
nm, E = 0.45|r| 7 = 10|r|, and 35 propagating channels. 



Appendix A: Nanoribbon in a 2DEG 
1. Unperturbed Green's function 

For a 2DEG confined in the x-dircction, creating a ribbon of width VF, the wave functions can be written as 

0„(rO=e*'=='^Xn(x), (Al) 

where r = (a;, y) and n is the mode number associated with the transverse eigenfunctions (assuming infinitely high 
confining walls at x = and a; ~ W) given by 



with the corresponding eigenenergies 



Xn{x) = y — sin(A:„a;), 



2[i 



(A2) 



(A3) 
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Here, fc„ = mr/W is the transverse momentum, ky the longitudinal momentum and fi the electron mass. Using these 
wave functions we may construct the free propagator, or Green's function, of an electron (having energy = E + irj, 
in the limit ?/ — > 0^) between the points f ' and r, in mode n, as 

a (rr'-F)- T ^MMlI 

^"^ ' ' ^"i-oo 2nE+-e^{ky) (A4) 
= Xn{x)xn{x')rn{y,y';E), 

where 

J-oo 27r E+ - ^{kl + kl) 
This integral can be evaluated using standard contour integration techniques)^ giving us that 



if E>E,„ 



rM;E)^{ (A6) 



and 



where En = {h^ /2p)kl 



^n{E)^^l'^\E-Enl (A7) 



2. Green's function, one impurity 

We introduce a single impurity modelled by an impurity potential with matrix elements Vnmi'r)- The perturbed 
propagator for an electron going from position f ' to r, while changing mode from m to n, can then be written using 
the Dyson equation as^ 

Gn7n{r,r';E) = gn{r,r';E)5nm 

+ X! 3n{r, E)Vni{fi)Gim{n,r'; E) 
J ^ , (A8) 

= gnir,r'; E)Snm + Gnm{r,r' ] E). 

Here we assume that the impurity is positioned at fi = {xi, yt) and that its potential is highly localized (delta-function 
shaped) so that all matrix elements but Kim(n) are zero. After introducing 

Turnin; E) = Vnm{ri) + ^ Vnl{fi)gi(fi, ff, E)Tim{fi; E) , (A9) 

; 

the scattering part of the Dyson equation can be rewritten on the T-matrix form 

Gnm{f,f';E)=gn{f,n;E)Tnrnin;E)g„^{n,f';E), (AlO) 

where 

Tn,n{n; E) = Vnmin) + ^ Vnl{?^)gl{n, n- E)Turr{n; E). (All) 

I 

Since the impurity is highly localized in position space, we may further assume that it will scatter equally between 
all different modes n and m and we have that Vnmi'f^i) = ^('^i) ~ 7 where 7 is the impurity strength. Using this 
assumption, it follows that T„m{f', E) = r(r; E) and we find that 



T{n-E) = Vlf,) + V{f,) 



7 



^gi{n,n]E) 



Tin;E) 



1 -iY.i9i{n,n;E) 
1 

ll-f + a,{E)+iap{Ey 



(A12) 
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where 

'^e/piE) -J^l^ ^ (A13) 

and e and p are the sets of all evanescent [E < Ei) and propagating (E > Ei) modes. Inserting the above expression 
for T(r,;; E) back into Eq. ()A10p allows us to solve for Gnmir, r'; E) and consequently for Gnmif, r '; E). 



by 



3. Fourier transformed density of states 

Once the perturbed propagator is known, the change in the local density of states (LDOS) due to scattering is given 



/5(r ; E) = Im Gnm [r, r; E) 

TT ^ — ' L 

nm 

IT f ^ 



(A14) 



where 



and 



h^J {I + a,{E)y + a^iE) k^{E)k^{E) ' 



-~fscXKniE),Km{E)) if n, TO € p, 

/es(K„(^),0)e--™(^)l?'-?'-l iinep,mee, 

fcs{0, K™(i;))e-«"(^)l^-^-l if n e e, TO e p, 

_-crp(i;)e-('^''('^)+«'"('^»l2'-2''l if n,TO e e, 



(A15) 
(A16) 

(A17) 



where 



and 



/cs(k1,K2) = (1 + 0-e)c0S + K2)\y ~ fji 

+ ap sin [{ki + K2)\y - yt\] , 

/sc(ki, K2) = (1 + <Je) sin [{ki + K2)\y - Vi 
- o-p cos + K2)\y - yi\] ■ 



(A18) 



(A19) 



When taking the Fourier transform of the scattering LDOS, we want to be able to resolve differences in x-momenta 
equal to or greater than -k/W (since this is the separation in kx, or fc„, between two adjacent subbands). This requires 
us to integrate over the interval \—W^ W] and we extend the function Pnmi^) such that it is even with respect to the 
origin. The Fourier transform is then defined as 



Afn7n{q;E) = K.„m{E)x 

00 

00 

dy 



n' — ~oo 



TT W 



w 



dx 



2tt 



e-^'-ypl^Xy^E), 



(A20) 



where the comb function fixes to multiples of tt/W. The x-part of the Fourier integral is 

dx _ 1 



(A21) 
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a) 



b) 




FIG. 12. (a) The geometry of the armchah nanoribbon. (b) The first Brillouin zone with three sets of Dirac cones. 



independent of if n and m are evanescent or propagating modes. 

The y-part will depend on mode types. We have already shown what happens when n,m £ p. In addition, if 

n,m Cz e we get that 



If n S p, m G e then 



7le-^''^ypy^iy;E) = [Sy,{nn{E) - qy , n.^{E)) + Sy,{nn{E) + qy , k^{E))\ 



2-n 



where 



27r 



SUa,b) 



{lh + ae{E))h + ap{E)a 



+ a'^ 

If n e e, TO G p we just need to interchange the n and to in the expression above. 



(A22) 



(A23) 



(A24) 



(A25) 



Appendix B: Armchair graphene nanoribbon 

In this appendix we first derive an analytic expression for the Green's function of an armchair nanoribbon with a 
single impurity. For the geometry, see Fig. Il^ a). We then derive the Fourier transformed density of states. 



1. Unperturbed Green's function 



The first Brillouin zone (IBZ) of graphene contains one pair of inequivalent Dirac cones. It is necessary, however, 
to include three pairs of cones [see Fig. HHb)] in order to incorporate all scattering events. The cones are located at 
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Kf = {±K^,0), = {±K^/2,Ky) and Kf = {±K^/2, -Ky), where K,^, = 47r/3a and Ky = 27r/3ao. The distance 
between two neighbouring atoms is denoted ao, while the lattice constant is denoted a = -v/Sao- This gives us three 
sets {d — 1, 2, 3) of wave function spinors^^ 

^dni^ = ilPf^}^ Xdnix), (Bl) 



where the longitudinal wave function components are 



= Ai^^i£|le^(^-+'^«)^ (B2) 



m^^{y) = ie'^^-'y+'^y^y , (B3) 

and the transverse wave function is 

Xdn{x) = 2Csin [{Kdx + kdn)x] . (B4) 

The integer number n labels the quantized transverse momentum kdn = rnr/W — Kdx in cone pair d. For each 



mode n, we have positive and negative energy subbands edn\{ky) — Xvf^k^ + {kdnY labeled by A = ±1. The Fermi 
velocity is Vf — 3ao|i|/2, where t is the nearest neighbor tight-binding hopping energy. The wave functions have been 
normalized through a normalization constant C = yJl/4W found from the condition dx\xdn{x)\'^ = 1/2. Thus, 
Xdn{x) = Xn{x) = ^JTJW s\n{mT/Wx). 

The free propagator for band n (in cone pair d) is computed as 

dky ^dn{^^l^{r') 



gdn(r,r ') = J2 

\ i_ 1 — 



27r E+ - tdn\{ky) 



where 



and 



xn{x)xn{x ) [rZHy,y';E) rf^^iy,y';E)J ' 



dL, e^fcyfe-?/') 



^dn [y,y,E)-T2tVfe I 2n {E+y-vj{kl + klJ- ^^^> 



After contour integration, we find the final form to be 



and 



where 



ri^/^^{y,y';E) = -^e^'<'^y(y-y'^ 

Vf 



zsgn (E) kdn 



± sgn {y - y') 



^issYv{E)K.a„(E)\y~y'\ 



(B5) 



(y, y,E)-ZEe / _ ^^^^^ ^ , (Bb) 



ri^/^^(2/, y'; E) = -zLI e^^^.C^^^ ) , (B8) 

Vf Kdn[-l^) 



(B9) 



'^dn = ^\iE/vfr-kjJ. (BIO) 

In the above formulas, we have assumed that n is a propagating mode (e.g. \E/vf \ > \kdn)- If mode n is evanescent 
{\E/vf\ < \kdn\), we have to modify the longitudinal momentum so that Kdn — ^ isgn(i?)Krf„. 
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2. Green's function, one impurity 



For the graphene armchair ribbon, we select an impurity fully localized on the A-sublattice, scattering equally 
between all modes n and m. The matrix elements of the impurity potential then is 



1 




(Bll) 



where 7 is the impurity strength. The T-matrix equation is written down in analogy to the 2DEG case, but for 
graphene it acquires a 2x2 matrix structure. For the potential in Eq. (jBll[) . we get 



. d I 

1 



T(f„r;) 



(B12) 



where 



and 



(B13) 



(B14) 



The letters e and p denotes sets of evanescent and propagating modes, respectively. The Dyson equation for the 
Green's function can now be written as 



where 

Gdncni{r,r'\E) = ^dn{r,n\E)T{n,ri;E)g,cm{n,r'\E) 
The scattering part Gdncm{T^,r''] E) takes the form 

1 



(B15) 
(B16) 



Gdncm{r,'r';E) = 



gf^{r,n;E)g^^{n,f';E) gl^{r,n;E)gi^{n,f';E) 

1/7 + <Te{E) + iap{E) \gf^{r, E)gi^{n,r'; E) gl^{r, r,; E)gf^^^^^^,. ^) 

Xn{x)Xn{x^)xMXm{x') (vj^ {y , y^,; E)T^^{y,,y' ; E) Tj^{y,y,-E)T^^{y,,y';E) 
1/7 + <ye{E) + i<y,{E) KT^^iv, y^■.E)T^^{y,, y'; E) rf„^(y, j/,; i?)rf^(y„ y'; E) 



(B17) 



For the computation of the local density of states, we need the two diagonal components. Their explicit forms 
(n, m (z p) are 



^dncmi^^ ^'\E) 



1 



\E\ 



lh + ae{E)+i(jp{E) \v) 



Xn{x)Xn{x{)X'm{Xi)Xva{x')x 



Jsgn{E){Ka„{E)\y-yi\+K,rr.iE)\y,-y'\) 
{E)Kcm{E) 



and 



l/j + ae{E)+iap{E) 



1 \^ 

- Xn{x)xn{x^)xMxm{x')e''''-'^y~y^^e'''^y'^y''y^: 



X gisgn(B)(K<j„(B)|y-yi|+Ke™(B)|j/i-y'|) 



isgn (E) kdn 

KdniE) 



sgn {y - Vi 



zsgn (E) kc 

>icm{E) 



(B18) 



sgn {yi - y') 



(B19) 
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3. Density of states 

The scattering correction to the local density of states can be computed separately for the two sublattices, and is 
given by 



TT 

dc nni 



(B20) 

-EE ^*(l(^)/5?i™(^)pi^^™^'(y; E), 



TT 

dc nm 



where Pnmi^) = Xn{x)Xm{x)- The A/B sublattice corrections are found by substituting Eq. (jBlSp and Eq. (|B19p 
respectively in Eq. (jB20l) . The results are very similar (the A correction being almost identical) to the 2DEG case, 
and for n, to € p we find that 

^lc™(^) = -(«//|i?|)'/Ci_(i?), (B22) 
Pdncmiv^ E) = <Jp{E)FSncm{y " y.; E) - (1/7 + a,{E))Fl,,^{y - y,-E), (B23) 



and 



Pdncmir^E) = Wp{E){-knk^ + n„{E)K„,{E)) 

+(f /7 + cTeiE))sgn (E) sgn {y - yi) (fc„K„(£;) + fc„K„(£;))} i^d"„„„(y - y^■,E) 
+ {ap{E)sgn (E) sgn {y - yi) {k {E) + k (E)) 

-(f /7 + (Je{E))sgn (E) sgn {y - yi) (fc„K„ {E) + fcjnKn {E))}F^^^„^{y-yf,E) 



(B24) 



where 



and 



FLcmiv, E) - COS [sgn{E){Kdn{E) + n,,r.{E))\y\ + {Kdy - K,y)y] (B25) 



FLcm^y. E) = sin [sgn{E){ndn{E) + K„r^{E))\y\ + [Kdy ~ K,y)y] . (B26) 

The Fourier transform of each component is carried out exactly as for the 2DEG, using Eq. (jA20[) . and the results 
for the AGNR are shown in Section ITVl 
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